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Magnetized Accretion Flows: Effects of Gas Pressure 



o 
o 

1—5 



Oh; 

o 



(N 
> 
O 
(N 

a^ 
rn 

in 
o 
a^ 
o 



M. Moscibrodzka^* and D. Proga^ 

^ Department of Physics, University of Illinois, 1002 West Green Street, Urbana, IL 61801, USA 

^ Department of Physics, University of Nevada, 4505 South Maryland Parkway, Las Vegas, NV 89154, USA 



Dates to be inserted 



ABSTRACT 

We study how axisymmetric magnetohydrodynamical (MHD) accretion flows depend 
on 7 adiabatic index in the polytropic equation of state. This work is an extension 
of Moscibrodzka & Proga, where we investigated the 7 dependence of 2-D Bondi- 
hke accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if 
simulations for various 7 can give us insights into to the problem of various modes 
of accretion observed in several types of accretion systems such as black hole binaries 
(BHB), active galactic nuclei (AGN), and gamma-ray bursts (GRBs). We find that for 
7 4/3, the fast rotating flow forms a thick torus that is supported by rotation and 
gas pressure. As shown before for 7 — 5/3, such a torus produces a strong, persistent 
bipolar outflow that can significantly reduce the polar funnel accretion of a slowly 
rotating flow. For low 7, close to 1, the torus is thin and is supported by rotation. The 
thin torus produces an unsteady outflow which is too weak to propagate throughout 
the polar funnel inflow. Compared to their HD counterparts, the MHD simulations 
show that the magnetized torus can produce an outflow and does not exhibit regular 
oscillations. Generally, our simulations demonstrate how the torus thickness affects 
the outflow production. They also support the notion that the geometrical thickness 
of the torus correlates with the power of the torus outflow. Our results, applied to 
observations, suggest that the torus ability to radiatively cool and become thin can 
correspond to a suppression of a jet as observed in the BHB during a transition from 
a hard/low to soft/high spectral state and a transition from a quiescent to hard/low 
state in AGN. 

Key vifords: accretion, accretion discs - black hole physics - X-ray: binaries - galax- 
ies: nuclei - methods: numerical -MHD 



1 INTRODUCTION 

Many astrophysical objects powered by accretion on to a 
compact objects like e.g. active galactic nuclei - AGN, black 
hole binaries - BHB, gamma-ray bursts - GRBs, show sig- 
natures of outflows, either in a form of wind or jet, or both. 
Large scale, coUimated jets are being detected fro m radio 
to 7- r ay w avelengths in AGN and microquasars (e.g. iKonigll 
I2OO6I . iMarkoff 2006 ). Recently an increasing radio emission 
from the core of some BHB was associated wit h small scale 
outflow/jet-like structures l|Fender et al.|[20o3 ). The obser- 
vations of BHB suggest that the power of the jet-like outflow 
may be related to the evolutionary stage of the object. 

The BHB exhibit several X-ray spectral states, which 
can change in time along a hysteresis loop on the luminosity- 
hardness diagrams. BHB (as well as galactic nuclei) spend 
most of their time in the quies cent /off mode du ring which 
their X-ray luminosity is low (|Gallo et al.ll2003l ). The dis- 
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covery of 17 out of 20 known BHB was possible only during 
an outburst (X-ray novae) accompanying the transition from 
the quiescent/off to a low/hard spectral state after which the 
X-ray brightness increases three or more orders in magni- 
tude (|Remillard fc McClintockll2006l '). In the low/hard spec- 
tral state the non-thermal power-law component dominates 
in the X-ray energy band. After the transition into the so- 
called high/soft state the spectrum is dominated by the ther- 
mal multi -blackbody emission co ming from a cold Keplerian 
disc (Sha kura fc SvunvaevlflQT^ ) . The source returns to the 
quiescent state through the soft-state, which closes the hys- 
teresis cycle. It is also well established that spectral states of 
BHB can change on r elatively short time-scales from days 
and months to years (i Remillard fc McCli ntoc^ \200dt ) and 
that the X-ray lumonisty of the BHB d oes not n ecessarily 
correlate with its spectral state ([Done et alll200'i1 ). 

The spectral state changes probably reflect the changes 
in the accretion rate/pattern and the formation/suppression 
of an outflow. For example, in the low/hard state, the BHB 
radio luminosity is strongly correlated with the X-ray lumi- 



2 M. Moscibrodzka and D. Proga 



nosity (Lr ~ L^ MGallo etlllbOOal ). This correlation does 
not hold after the transition to the hi Rh/soft state, when 
the radio em ission becomes suppressed (jCorbel et al.|[200ll . 
iGallo et alj [20031 . The X-ray variability and radio activ- 
ity strongly support a model, where in the low/hard state, 
a standard thin accretion disc is truncated and does not 
reach the black hole horizon. During low/hard state, the 
inner disc that used to be during the high/soft state is re- 
placed by a hot thick flow, which probably acts as a lunch- 
ing site for a steady j et with the Lorentz factor, F, typi- 
cally less than 2. (e.g. 'Heinz fc Sunyaev 'ioO^; jMeied l2005l : 
iRemillard fc McChntoc k 2006; Fonder ct al. 200l and refer- 
ences therein). 

The transition from low/hard to high/soft spectral state 
occurs via an intermediate state, the so-called very high 
state. The very high state is characterized by strong black 
body and substantial power-law emission. In the very high 
state, the radio emission exhibits flaring. The flaring can be 
due to a fast non-stati onary optically thin outflow with high 
Lorentz factor, F > 2 (|Fender et al.ll2004 ). 

A radio- X-ray correlation, similar to that seen in BHB, 
is also expected to work in much larger scale systems like 
AGN. Because of much longer evolutionary time-scales in 
AGN, in comparison to BHB, the problem may be studied 
statistically only. The expectation of the correlation is sup- 
ported by the fact that the radio-loud quasars have higher 
X-ray fluxes than radio-quiet quasars (^Zdziarski et al.iil995i) . 
In this scenario, the AGN counterparts of the 'quiescent' 
state are Low Luminosity AGN (LL AGN), because they 
lack the 'big blue bump' in their spectra, are weak in X- 
rays and are loud in radio in comparison to AGN (e.g. IHoI 
I2OO3I : iKawakatu et al.|[2009h . Our own galactic center - Sgr 
A* - is considered as an extreme example of LL AGN. The 
models of the quiescent/off phase of AGN and BHB, are 
usually based on the two assumptions: either most of the ac- 
cretion energy is transferred into the outflow kinetic energy 
which results in a reduction of the mass-accretion rate. Ma, 
and l ower luminosity , or the accretion is highly inefficient, o r 
both (|lchimarulll977l . lNaravan et al.lll998llYuan et al.ll2003h . 

From the theoretical point of view, it is still unclear 
what induces the accretion phase transition in BHB, AGN, 
and GRBs. Many authors considered 1-D, vertically inte- 
grated, models where the cold disc evaporates, t ransitions 
into a hot phase and then rc-conde nsates (e.g 
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Mever et alli2000. Roiariska & Czernvll200d . 
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Mever et al.ll2007l ). In 1-D models, the mass-accretion 



rate is the key parameter responsible for the disc behaviour. 
These models include such physical processes as mass trans- 
fer between hot and cold phase, thermal conduction, and ra- 
diative transport. However, the 1-D models do not include 
an outflow formation and do not yet attempt to reproduce 
observed correlation between X-rays and radio luminosity 
probably c aused by outflo w, during low/hard state. As un- 
derlined in lAbramowicz et aD (j2000h . at least a 2-D treat- 
ment of the problem is needed to correctly model the accre- 
tion state transition process. 

In this paper, we present results of 2-D axisymmetric 
magnetohydrodynamical (MHD) simulations of a slowly ro- 
tating accretion flow on to a black hole (BH). We investigate 
the time evolution and properties of inflows/outflows for var- 
ious values of an adiabatic index, 7 . Our w ork is an extension 
of work of iMoscibrodzka fc Progal ^200^ (hereafter MP08) 



who considered models for var ious 7 but in the hydro- 
dynamical (HD) limit (see also IProga fc BegelmanI l2003al . 
he reafter PB03a, for HD sim ulations with 7 — 5/3, work 
of lProga fc BegelmanI [2003bl . hereafter PB03b, who consid- 
ered MHD models only for 7 — 5/3, and finally work of 
iProga et al.l (|2003l ) who considered a MHD coUapsar model 
for GRBs where the inner accretion flow is radiation pres- 
sure dominated and 7 is effectively « 4/3). 

In MP08, we showed that in the HD limit, many proper- 
ties of the accretion ffow strongly depend on 7. Most relevant 
to this work is that the thickness of the torus decreases with 
decreasing 7 whereas Afa in units of the Bondi rate, Mb 
increases with decreasing 7. Here, we want to check how 7 
affects magnetized flows. For example, is the efficiency of 
producing outffows sensitive to the thickness of the torus? 
This question is related to the observations of BHBs and 
AGN. One could argue that magnetically driven outffows 
should be the strongest from geometrically thick discs/torii 
with significant pressure support and the weakest, if any, 
from geometrical thin discs w ith little pressure support (e.g. 
iLivio et al.lll999l : iMeierlbOOll ). 

To our best knowledge, there are no direct simulations 
that would confirm the above expectation mainly because 
previous simulations considered thick, pressure/rotation 
supported torii. What is clear is that most of these pre- 
vious global MHD simulations of acc r etion ffows show 
outffows (e.g. 'Uchida fc Shibata' 'l985'; 'Stone & Norman] 
1994; Stone & Pringlc 2001; Hawlcv fc Balbus 2002; PB03b; 
Proga et a l. 2003 :D e Villiers et al. 2005; Mi zuno et al.l|2003 : 
Kato et all ,2004: .McKinnev fc Gammie .20041 : IProgal lioosl . 
Fragile fc Meied l2009l ). However, it is unclear under what 



conditions (e.g., for what torus thickness) the torus stops to 
produce an outffow. 

Our very simplistic approach to control the torus thick- 
ness is dictated by a complexity of the relevant physics. 
A full treatment of an accretion disc, especially geomet- 
rically thin one, requires modeling energy dissipation, ra- 
diative heating and cooling, and radiative transfer. In our 
approach, for 7= 5/3 and 4/3 the torus will be hot and sup- 
ported against gravity by gas pressure and rotation. As such 
it will be geometrically thick. On the other hand, for 7 ~ 1 
the torus will be supported only by rotation and therefore 
it will be Keplerian and geometrically thin (or at least con- 
siderably thinner than for 7 = 5/3). Therefore, our study 
of effects of 7 will allow us to gain some insights into the 
problem of the outffow formation for various physical condi- 
tions and different phases of activity like e.g. soft and hard 
spectral states in BHB and AGN activity modes. 

Here we will not be able to immediately compare our re- 
su lts with observations in an y great detail (as it was possible 
in IMoscibrodzka et al. I I2OO7I where we computed and com- 
pared radiative properties predicted by PB03b's MHD simu- 
lations with observations of Sgr A*). However, our approach 
enable us to produce results that can be easily and consis- 
tently compared with a great body of data from complemen- 
tary numerical simulations that were preformed within one 
framew ork and with the s ame numeri cal tools (i.e.,PB03a; 
PB03b: iProga et al.ll200i |Progall2005l , MP08). 

The article is organized as follows. In §[5] we describe our 
method and initial conditions used to set our simulation. In 
§[3]we present the results of the numerical MHD simulations. 
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The results are discussed in the contest of previous work in 

§H1 



2 METHOD 
2.1 Equations 

To calculate the structure and evolution of an accreting flow, 
we solve the equations of MHD 



Dp 
Dt 



+ pV -v = 0, 



(1) 



P^ = -VP + pV$ + ^(VxB)xi?, (2) 



D 

'di y-p 



= -PV ■ u, 



dB 

dt 



= Vx{vxB) 



(3) 



(4) 



where p, P, v, e, $ and B is the mass density, gas pres- 
sure, fluid velocity vector, internal energy density, gravita- 
tional potential, and magnetic field vector, respectively. La- 
grangian derivative is defined as: = - ^ . We relate 

P and e through the equation of state P = (7 — l)e, where 
7 is an adiabatic index. 

We perform simulatio ns using the pseudo-New tonian 
potential $ introduced by iPaczvnski fc Wiital (|l980l ) (here 
after PW) 



GM 



(5) 



where Rs is the Schwarzschild radius. This potential ap- 
proximates general relativistic effects in the inner regions, 
for a non-rotating BH. In particular, the PW potential re- 
produces the last stable circular orbit at r = 3Rs as well as 
the marginally bound orbit at r = 2i?s- 



2.2 Initial conditions and boundary conditions 

Generally, the setup and numerical methods used in our sim- 
ulations are as in PB03b. Here, we only briefly summarize 
the initial and boundary conditions and spell out the differ- 
ences between our and PB03b simulations. 

Our calculations are performed in spherical polar coor- 
dinates (r, 6, (f)). We assume axial symmetry about the rota- 
tional axis of the accretion flow {6 — 0° and 180°). For the 
initial conditions of the fluid variables we follow PB03b and 
adopt a Bondi accretion flow. In particular, we adopt Vr and 
p computed using the Bernoulli function and mass accretion 
rate for spherically symmetric Bondi accretion with the PW 
potential. We set poo ~ 2.2 x 10"'^'^ g/cm^ and mass of the 
central black hole Mbh = 3.7 x W^'Mq which make our re- 
sults directly applicable to the galactic center (e.g. it will 
be useful in our future calc ulations of synthetic spectra as 
in iMoscibrodzka et al.l [20071 ). However, our results could be 
rescaled and applied to other systems. We specify the sound 
speed at infinity, cSoa through a normalized Schwarzschild 
radius Rg = Rs/Rb {R's = 2cs^/c^), where Rb is the Bondi 
radius. The normalized BH radius characterizes the gas tem- 
perature. We consider models with R'^ = 10~^. 



We specify the initial conditions by adopting a non- 
zero specific angular momentum I at the outer boundary. 
We consider a case in which the angular momentum at the 
outer radius To depends on the polar angle via 



l{ro,0) = loil-\ cos e\). 



(6) 



This formula was used by PB03a and MP08. We note that 
PB03b use a step function so that this is one of the dif- 
ferences between the setup of our and PB03b's simulations. 
However, as discussed by PB03b and shown below the de- 
tails of the I distribution are not very important. We express 
the angular momentum on the equator as 



^0 



-.RbCSo 



(7) 



where R'q is the 'circularization radius' on the equator in 
units of _Rb for the Newtonian potential (i.e., GM/r^ = v^/r 
at r = R'cRb)- 

To generate the magnetic field, we use vector potential, 
A, i.e., i? = V X A. The initial magnetic field configura- 
tion is a pure vertical magnetic field. The zero divergence 
vertical field is defined by the potential A—{Ar = 0, Ae = 0, 
A^ — Ar sin 6). We scale the magnitude magnetic field by 
plasma parameter /3in = SnP/B^ at the inner boundary, ri. 



A^ \/CP~AM, 



(8) 



where P,-. is the initial pressure at the inner radius. 

Our standard computational domain is defined to oc- 
cupy the radial range n = 1.5 -Rs ^ J" ^ ''o ~ 1.2 Rb 
and the angular range 0° ^ 6 ^ 180°. The r — 9 domain 
is discretized into zones with 140 zones in the r direction 
and 100 zones in the 6 direction. We fix zone size ratios, 
drk+i/dr-k = 1.05, and dOi/dei+i = 1.0 for 0° < 6* < 180°. 
With such defined grid we are able to resolve the fastest 
growing mode of the magnetorotaitonal instability (MRI) 
only in a case of thick torus (7=5/3). See 13.61 

The boundary conditions are specified as follows. At the 
poles, (i.e., ^ = 0° and 180°), we apply an axis-of-symmetry 
boundary condition. At both the inner and outer radial 
boundaries, we apply an outflow boundary condition for all 
dynamical variables. As in PB03b, to represent steady con- 
ditions at the outer radial boundary, during the evolution of 
each model we continue to apply the constraints that in the 
last zone in the radial direction, ve = 0, = lof{9)/ro sinO 
(where f{9) = 1 — | cos ^|), and the density is set to poo at all 
times. Note that we allow Vr to float. For the magnetic fleld, 
at the outer boundary, we ap ply an outflow conditio n. At 
the inner boundary, we follow IStone fc Pringld l|200lD and 
use negative stress condition (i.e. we enforce BrB^ ^ at 
ri). To reduce the problem of high Alfvenic velocities in re- 
gions of low density, we set a lower limit to the density on 
the grid as pmin = p{r\lr)^^^ and enforce it at all times in 
the simulations. 

To solve Eqs. dlT l- Q we use the zeus-2D code described 
bv IStone fc NormanI l|l992al lbh. modified to implement the 
PW potential. 



3 RESULTS OF THE SIMULATIONS 

The parameters used for the simulations are summarized in 
Table [1] The Table columns (l)-(8) show respectively the 
number of a given run; the BH radius and the Bondi radius 
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ratio (-Rs)i ^^'^ circularization radius compared to the Bondi 
radius (-Rc); the adiabatic index (7); the initial magnetic 
plasma parameter (/3in) at To; the final time (if, in) at which 
we stopped each simulation in units of the Keplerian orbital 
time at rt (idyn ~ 595s for Mbh = 3.7 x IO^Mq); the same 
final time if, out but in units of the Keplerian orbital time 
at To, and finally the time-averaged mass accretion rate on 
to the BH in units of the Bondi accretion rate {Ma,/ Mb) 
over 1000 idyn at the end of each run. The shorter duration 
of our simulations (compared to the duration of PB03b's 
simulations) is caused by the high density gradients which 
appear in our runs for low 7, and which lead to significantly 
smaller time steps. 

3.1 Velocity field and sonic surfaces evolution 

In our simulations, the evolution of the poloidal velocity field 
initially proceeds in a similar manner as in the correspond- 
ing HD simulations (see e.g. panels 3 and 4 in fig. 2 in PM08 
for each 7). Namely, after a short episode of the radial in- 
flow, the matter with relatively high- / forms a torus in the 
equatorial plane around the BH (high-/ is defined here as 
I > Icr and low-/ as / < Icr, where Icr = 2Rsc is a critical 
specific angular momentum). For 7=5/3, 4/3, and 1.2, ini- 
tially when the torus forms, a shock wave is produced. The 
shock wave forms because some of the gas forming a torus 
turns around and moves outward along the equator due to 
the centrifugal force and gas pressure or tries to move closer 
to the poles and accrete. For 7 =1.01, the gas is the most 
gravitationally bound compared to other 7's. As such it ac- 
cumulates near BH and does not fiow out along the equator 
contrary to the higher 7 cases. 

In its early phase of the evolution, the torus is not able 
to accrete on to a BH, independently of 7, because the MHD 
turbulence has not yet fully developed to transoport the 
angular momentum. The matter in the polar regions does 
not feel any centrifugal barier, thus it flows radially toward 
the BH. In our MHD simulations and the HD simulations 
studied in PM08, the sonic surfaces topology evolves in a 
similar manner. It changes from a circular shape into two 
lobes where the lobes are elongated along the poles. As we 
reported in PM08, the evolution of the sonic surface is the 
slowest for 7=1.01, due to a lack on an equatorial outflow. 

After a short time the magnetic field effects start to 
dominate the accretion fiow evolution. The torus starts to 
accrete toward the BH. The biggest and expected differ- 
ence between HD and MHD simulations is the development 
of bipolar outflows in MHD simulations regardless of 7. 
PB03b's simulations well illustrate the general evolution of 
the magnetic flelds. The polar outflows form at small radii 
(r < 20-Rs) close to the poles and expand to the outer 
boundary of the computational domain along the symme- 
try axis. We flnd that the time at which the outflows form 
depends on 7: the lower 7, the later the outflow forms. Al- 
though initially HD and MHD flows evolve similarly, at the 
end of the simulations, the geometry and character of the 
accretion flows in the HD and MHD simulations, are sig- 
niflcantly different. At the end of the simulation, the MHD 
flows are composed of a bipolar outflow propagating in the 
polar regions, and the inflow near the equator. The equa- 
torial region is turbulent and strongly asymmetric for all 
7's. 



3.2 Mass accretion rate evolution 

Fig. [T] shows the Ms./Mb as a function of time, for each 
studied 7 as listed in Table [T] The time is given in units 
of the orbital time at the inner boundary of the computa- 
tional domain. A pair arrows in each panel mark the mo- 
ments at which the polar outflow forms, in the Northern 
(N) and Southern (S) directions. The moment of the forma- 
tion of the large scale polar outflow is defined as a moment 
when the matter stops accreting through the polar regions 
and the radial velocity sign changes sign from negative to 
positive. We determine this by inspecting the velocity maps 
and checking that the outfiow velocity becomes larger than 
escape velocity at the outer boundary. 

The upper-left panel in Fig.[T]shows Ma/Ms for run 1. 
Initially, Ma /Mb decreases from 1 to 0.3 and then remains 
constant until t ~ 2500. The value for Ms./ Mb of 0.3 is con- 
sistent with that obtained in the HD counterpart (see MP08 
and PB03b), and it corresponds to the matter accreting in 
the polar funels in the initial HD phase of the simulation. 
However, later around t= 3700, Ala./ Mb suddenly decreases 
to a minimum value (below 0.001), and oscillates between 
the minimum value and 0.03. We note that the moment of 
the S outfiow formation does not correspond exactly to the 
moment of the Ma/ Mb reduction to its minimum value. In 
general (the details will be discussed in the next sections) 
the mass accretion rate reduction is caused by the men- 
tioned bipolar outflow, which pushes the radially infalling 
material away from the polar funels. Such a vacuuming of 
the polar regions makes the torus the only source of gas 
available for accretion. During later phases of the evolution 
(i.e., t > 4000), Ma/Me shows narrow spikes and dips simi- 
lar to those found by PB03b, corresponding to the episodes 
of enhanced and suppressed torus accretion, respectively. 

The upper-right panel in Fig.[T]shows the Ma/ Mb evo- 
lution for 7=4/3 (run 2). At the beginning but after the 
torus formation. Ma/ Mb ~ 0.2, which is consistent with 
mass accretion rate found previously in the HD simulations. 
At this stage, the gas is being accreted from the polar re- 
gions. Later when an outflow forms at i « 4000, Ma /Mb 
suddenly decreases. For t > 4000, the Ma/Ms evolution is 
characterized by narrow spikes, dips (similarly to run 1), 
and additionally broader features with higher amplitudes. 
The mass accretion rate in the spikes can reach a level of 
about 0.02 and corresponds to the torus accretion. During 
the dips, the Ma/Me has a minimum value below 0.001, 
when the torus stops accreting. The rates seen in the two 
broad features (around t=5500 and t=7000) reach the level 
of ~ 0.1. At even later times {t > 9000), Ma/ Mb is almost 
constant around 0.01. The broad features have been found 
before by PB03b, and correspond to the non-rotating gas ac- 
cretion that reenter the BH vicinity in the equatorial plane. 
Near the end of the simulation, the minima (dips) reappear. 

For run 3 (the lower- left panel in Fig. [T|, the evolution 
of the mass accretion rate differs signiflcantly from runs 1 
and 2. In particular, there are no sudden drops in Ma/Ms- 
Instead, Ma/Me evolves quite smoothly, despite the forma- 
tion of bipolar outflows. The Ma/Me curve is variable, and 
it shows two significant minima around t=7000 and at the 
end of the simulation. The mean level of Ma/ Mb is around 
0.1, which is rather high compared to those of runs 1 and 2. 

The lower-right panel in Fig.[l] presents Ma/ Mb for the 
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Table 1. Summary of parameter survey of the MHD simulations. 



run 








7 






tf.out 


A/a/MB 


1 


10-3 


5 X 10- 


-2 


5/3 


109 


1.22 X lO"' 


0.55 


0.03 


2 


10-3 


5 X 10- 


-2 


4/3 


109 


1.2 X 10* 


0.54 


0.01 


3 


10-3 


5 X 10- 


-2 


1.2 


109 


8.76 X 103 


0.4 


0.1 


4 


10-3 


5 X 10- 


-2 


1.01 


109 


1.36 X 10"' 


0.62 


0.35 



model with 7=1.01 (run 4). Here, Ma/Ms shows a strong 
time variability although the time averaged rate remains 
high compared to other models. There is one significant 
minimum in Mb,/ Ms at the final stages of the simulation. 
The panel shows many periods of increased accretion, during 
which Ma /Mb is higher than 1! 

Interestingly, for 7=1.2 and 1.01, in spite of the presence 
of the bipolar outflows, Ma/Ma does not decrease as sharply 
as for 7=5/3 and 4/3. It also does not show the characteristic 
features such as spikes and dips which were also found by 
PB03b. The different character of the mass accretion rate 
evolution in runs: 1, 2 and 3, 4, reflects two distinct outflow 
formation mechanisms. 

3.3 Inflow/outflow rates 

To investigate the mechanisms responsible for the polar out- 
flow formation, we examine inflow/outflow/total fluxes of 
various physical quantities as functions of the distance from 
BH. The fluxes are averaged over the whole 9 range and over 
1000 tdyn before the end of each run. 

We compute the fluxes of: mass, internal energy, kinetic 
energy, magnetic energy and total energy. The mass flux at 
a given radius, r is given by: 

M ^r^ d) pvrdO, (9) 

where dQ = sm{9)d9d(l>, p is the density, and Vr is the radial 
velocity. To calculate the mass inflow rate Min at a given ra- 
dius, we include only contributions with Vr < 0, whereas to 
calculate the mass outflow rate Mout, we include only con- 
tributions with Vr > 0. To compute the net (total) mass flow 
flux Mtot, we include all contributions regardless of the 
sign ( Mtot = Min + Mout)- To compute the internal, kinetic, 
magnetic and total energy fluxes, we use the same approach 
as for the mass fluxes with small modification. Namely, p 
in Eq. [9] is substituted by e, Skin = p{vr + Vg + ^0)/2, 

S'poynting = (v X B)/Vr and Etot = i^kin + C + Spoynting 

for the internal, kinetic, magnetic and total energy, respec- 
tively. We normalize all energy fluxes with Mbc^. 

In Fig. [2] and Fig. [3] we only show the total mass flux 
and total energy fluxes, respectively. We present it as a 
log{l + Mtot) function so that the positive values will in- 
dicate the domination of an outflow and the negative - an 
inflow. The solid, dotted, dashed and dot-dashed lines cor- 
respond to runs 1, 2, 3, and 4, respectively (as marked in 
the panels). Below, we will also the main conclusion which 
will be based our analysis of all studied fluxes not just those 
presented in the flgures. 

The total mass flux (Fig. [2]) as a function of radius is not 



constant. This indicates that overall the flows did not reach 
a steady state. However, in run 1 and 2, the total mass flux 
is constant in the inner parts of the computational domain 
only up to few Rs- For models 3 and 4, the flux is nearly 
constant over a wider radial range, i.e., up to r ~ 100i?s- 

Our analysis of the various fluxes showed that in run 
1, the total internal energy flux changes sign from negative 
(domination of the inflow flux) to positive (domination of 
the outflow flux) very close to the BH (r ~ 10i?s). The 
similar behaviour characterizes run 2 and 3. Run 4 differs 
from other cases: its total internal energy flux is negative for 
all radii. 

There are a few general trends in the behaviour of all 
fluxes in all runs. The strongest contribution to the total 
energy inflow flux (upper panel in Fig. O in the outer re- 
gion of the computational domain (which we define here as 
the region for r ^ 10J?s), is from the internal energy, re- 
gardless of the value of 7, whereas the total energy outflow 
flux (middle panel in Fig. |3]) is dominated by kinetic energy 
(with one exception for 7 = 1.01 for which the internal en- 
ergy outflow flux is even greater than kinetic one but orfly 
up to r ~ 100i?s)- 

In the inner part of the flows (for r ^ 10i?s), where the 
most of the outflow is produced, the total energy infiow fiux 
is dominated by the kinetic energy for all runs. 

The total energy outflow flux in the close vicinity to 
the BH (i.e. for r ^ 10i?s) is dominated by magnetic energy 
in runs 1, 2, and 3. In run 4 (7 = 1.01), the total energy 
outflow fiux is dominated by the kinetic and internal energy 
- the contributions from both components are comparable. 
This suggests that the centrifugal forces may be responsible 
for an outfiow in run 4, rather than magnetic forces (runs 
1-3). We are elaborating this point below. 

3.4 Formation of a stationary outflow 

In this section, we investigate how the large scale bipolar 
outfiow is formed. W e will do it in the conte xt of previou s 
work, namely PB03b. IProga et al.1 l|2003D and lProgal l|2005l ). 
We will concentrate on identifying the fiow components and 
checking if they are by nature the same ones as those found 
in previous simulations. Therefore we first briefly summarize 
the results from PB03b. 

As mentioned in the introduction, PB03b considered 
models similar to ours but only for 7 = 5/3. They found that 
the nature of a weakly magnetized, slightly rotating accre- 
tion fiow is controlled by magnetic fields when I 2Rsc. As 
in the HD inviscid case, studied in PB03a (see also MP08), 
the magnetized material with I <; 2i?sc forms an equato- 
rial torus. However, in the MHD case, the torus accretes 
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on to BH because of m aKnetocentrifuRal instability (MRI; 
iBalbus fc HawlevI IToQll ) and magnetic braking. The mag- 
netic field in the torus amplifies and produces both a mag- 
netic corona and an outflow. The latter two can be strong 
enough to prevent accretion of the low-Z material through 
the polar regions, which was the source of accretion in the 
HD inviscid case. 

Very similar time evolution and flow behaviour was 
found in simulati ons of the MHD coUapsar model for 
gamma-ray bursts ijProea et al.ll2003l 'l which were similar to 
those of PB03b's simulations, except for much more sophisti- 
cated micro-physics. In particular, Proga et al.'s considered 
an equation of state that included contributions from gas, 
photons and electrons. The most relevant aspect of Proga et 
al.'s simulations for us here is that the inner torus is radia- 
tion pressure dominated and the effective adiabatic index of 
the inner torus is ~ 4 /3. A general conclusion from PBOSb 
and IProea et al ] |2003l is that for 7^4/3 the torus forma- 
tion and evolution, including the development of a corona 
and outflow are qualitatively similar. In particular, the torus 
outflow is strong enough to suppress polar funnel accretion 
and to propagate out to a large distance from the BH. This 
result is perhaps unsurprising given that the torus is geo- 
metrical thick for 7 ^ 4/3. 

The previous simulations showed that slightly rotating 
accretion flows are very variable and complex. In particular, 
the inner radius of the torus and the strength of the out- 
flow change with time. For these re asons, only a l ater, much 
more detailed analysis of PB03b and lProga et al.l (|2003l ) sim- 
ulations showed an additional flow component: a very col- 
limated outflow driven from an infalling gas with / ^ 2Rsc 
(|Progal [2005! ). To articulate the basic p hysics t hat o ccurs 
in production of this collimated outflow, IProeal (|2005l l per- 
formed also some new simulations. He considered gas with 
I < 2Rsc so that after an initial transient, the flow in the HD 
case accretes directly on to a BH without forming a torus. 
However, in the MHD case, even with a very weak initial 
magnetic field, the fiow settles into a configuration with four 
components: (1) an equatorial infiow, (2) a bipolar outfiow, 
(3) polar funnel outflow, and (4) polar funnel inflow. The 
second flow component of the MHD flow represents a sim- 
ple yet robust example of a well-organized inflow/outflow 
solution to the problem of MHD jet formation. 

In each of our present simulation, the accretion flow 
forms a torus and outflow although some low-/ {I < 2Rsc) 
inflow may occur for most of the time. We conflrm that an 
outflow may form in two manners. First type of an outfiow is 
driven by ce ntrifug a l and magnetic forces as it was already 
described in IProga (|2005l ). The second type of outfiow is a 
wind driven from a torus by the toroidal magnetic field gra- 
dient as described in PB03b (2003b and references therein). 
In the first case, the outfiow is collimated and removes only 
a small fraction of the low-Z matter from the polar funnel 
and does not affect the torus evolution, i.e. the torus is sep- 
arated from the collimated outfiow by an accreting stream 
of low-/ matter. In the second case, the outfiow is a broad 
wind from the torus. This wind can remove all low-/ mat- 
ter from the polar funnel (by pushing it sideways) . The two 
types of outfiows contribute the time variability of Ma/Me 
described in the previous section and of the fiow structure 
detailed here. 

In Figs.|4]and[5l we show the two-dimensional structures 



of the density and velocity fields along with the contours of 
the angular momentum of the gas. We also show the maps of 
the magnetic field parameter f3 — SirP/B^, including total 
magnetic filed {B^ = + Bg + -B|), at the times of the 
formation of an outflow on N and S side. 

For 7=5/3 (run 1), the high-/ matter (/ > 2Rsc, dotted 
lines on the angular momentum contours) does not reach 
the BH and a torus does not form yet ( here, we deflne the 
torus as a high- / matter accreting in the equatorial region). 
Instead in the inner fiow, the low-/ gas and the high-/ gas 
are mixed. Such a state of accretion was not found in the 
previous simulations because here we assume a continuous 
distribution of / with ^ at ro, whereas in the previous simu- 
lations of PB03b a step-function was assumed. 

In run 1, the initial outfiow is driven by the centrifu- 
gal and magnetic forces along the rotation axis as found by 
[Proga (2005) (but not in PB03b). As we mentioned in the 
previous subsection Ms./ Mb does not decrease significantly 
during the formation of the initial outfiow (Fig. [1] left-upper 
panel) , but it does decrease slightly some time after the out- 
fiow is formed. 

The polar outfiow is highly collimated by the infalling 
ambient gas. The outfiow that forms on the opposite side of 
the torus later during the simulation is very similar to that 
just described. At the moment of the outfiow formation, the 
torus is not developed yet. The outfiow is marked by the 
solid contours in the angular momentum panel in Fig. [3] 
(upper middle panel). As one can see this outfiow is made 
of the slowly rotating matter near the symmetry axis. 

In the later phases of the simulations (Fig. [5} , after the 
bipolar outfiow expansion, the high-/ gas forms a well de- 
veloped torus that extends inward to the last stable orbit 
and accretes on the BH. Subsequently, smaller scale (within 
20 -Rs) magnetic torus corona and outflow develop. In this 
phase, the torus outflow starts to push aside the low-/ mat- 
ter accreting from the polar regions. For the time up to t« 
6000, the high- and the low-/ gases mix with each other, 
near the BH. Later, when the low-/ matter is signiflcantly 
pushed aside by the fully developed magnetic torus corona 
and outflow, the accretion flow becomes very similar to that 
in PB03b (see also Sect. l33)) . 

For 7=4/3 (run 2, second row of panels in Figs.|4]and[5]), 
the outflow is formed in the similar way as found by PBOSb 
i.e. it is a collimated wind from a torus magnetic corona. 
In this case, the torus is more bound than in the case of 
7=5/3. Therefore the high-/ matter does not mix with the 
low-/ matter near the equator close to the BH (the upper 
panels of fig. 8 in MP08 for the HD run 1). The magnetic 
torus corona and outflow form relatively fast. The corona in 
this case extends up to r ~ 70-Rs- Because the torus is geo- 
metrically thick, and the low-/ matter is pushed aside very 
efficiently from the polar region, the polar infiow is com- 
pletely suppressed ( see Fig. [S] middle panel corresponding 
to 7 = 4/3). The moment of the sudden decrease of Ms,/Mb 
(Fig. [T| corresponds to the time of the outfiow formation. 

Next, we describe what happens for 7=1.01 ( run 4, 
bottom panels in Figs. |3]and[Sl). Already in the HD simu- 
lations (runs G, H, I, J in fig. 8 in MP08), one may notice 
that because of much lower gas pressure, the torus is ge- 
ometrically thinner compared to the cases with higher 7. 
Here, a strongly turbulent magnetic corona forms above the 
disc. The corona is geometrically thick, but it does not ex- 
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pand into the polar regions as in the previous cases. As in 
run 1, initially the outflow is caused by the cent rifugal and 
magnetic forces on both sides (as in IProgal l2005h . However, 
this colhmated outflow does not completely prevent the \ow-l 
matter inflow toward the BH from the polar regions. There- 
fore, Mb./ Ms remains always high (close the Bondi rate). 
In the lower panels in Figs. |4] and [5] the dotted contours 
in the angular momentum panels indicate that the high-/ 
matter is separated from the outflow by the low- / accreting 
matter marked as the solid contours. The low-/? corona is 
separated from the low-/3 collimated outflow with the high- 
P region of the accreting matter. As we mentioned in the 
previous section, the moment of the outflow formation can- 
not be clearly identifled in the mass accretion rate curve 
(Fig. [T}. The persistently accreting low-/ matter separates 
the polar outflow and magnetic corona for most of the sim- 
ulation time. In this case, the accretion flow consists of four 
parts (from the equator to the pole): (1) geometrically thin 
torus in the equatorial plane, (2) highly turbulent magnetic 
corona, (3) low-/ matter inflow through polar funnels, and 
(4) polar outflow. 

In the intermediate case of 7=1.2 (run 3, third row of 
panels in Figs. |4] and [5}, although the magnetic corona forms, 
the polar outflow forms in both manners as described for 
runs 1 and 2. In Fig. |4l one notices the corona formation on 
the N and S sides. The corona extends beyond r = 100-Rs, 
which is clearly seen in Fig. 2] First, on the S pole, the 
expanding corona causes the outflow. The outflow on the 
opposite pole, is c aused by the magneto-centrifugal mech- 
anism described in iProgal (|2005l ). The mass accretion rate, 
just as in the case of 7=1.01, does not experience a sudden 
reduction during its evolution because the outflow is weak 
so that the low-/ matter persistently accretes and mainly 
contributes to the total mass accretion rate. Nevertheless, 
the mass accretion rate for 7=1.2 is slightly lower than that 
for 7=1.01. 

In our models, the coUimation of the jet-like outflow 
depends on the 7 index and it is related to the mechanism 
of its formation rather than the torus geometrical thickness. 
For 7 = 5/3,4/3 and 7 = 1.2 (S side), the jet-like outflow 
at the flnal time of the simulation is not well collimated, 
because it is driven from the torus corona. To investigate the 
problem of coUimation, one should perform the simulation 
over a longer time period to allow the flow to reach the 
quasistationary state at larger radii. For 7 = 1.01 and 7 = 
1.2 (N side) the outflow is more collimated in comparison to 
the coronal jet-like outflow because it is formed by the low-/ 
matter and collimated by the accreting low-/ matter. 

3.5 Episodic coronal outflows for low 7 

For 7=1.01 (as well as for 7=1.2), M^jM-g, remain high for 
most of the simulations time. However, at the late phases of 
the simulations, we observe a few signiflcant minima in the 
mass accretion rate curves (lower panels in Fig[T}. These 
minima are caused by small scale coronal outflows which 
temporary suppresses accretion of the low-/ matter in the 
polar funel. 

Fig. |6] presents the sequence of maps of the logarithmic 
density, angular momentum contours, and plasma parameter 
/3, before, during, and after an episode of an outflow in run 
4. To better show the structure of the inner flow, we present 



only a half of the computational domain above the equator 
within 500 i?s. Before the outflow, the magnetic torus corona 
extents to 100 iis (the two top rows of panels). During an 
outflow, the corona expands up to 300 iis (the third and 
fourth rows of panels). Later, the matter flows back toward 
the BH. Before and after the outflow, the torus corona and 
polar outflow are separated by a stream of the low- / matter. 
After the outburst, Ms./Mb increases back to its mean high 
level (the bottom row of panels) . The whole cycle lasts about 

3500 tdyn. 

3.6 The structure of the torii 

In this section, we return to our investigation of the time 
evolution of Ma/Ms for various 7. Here, we identify accre- 
tion states that correspond to various episodes of torus ac- 
cretion. Figs.[7]shows the segments of the M^./ Ms evolution 
for runs 1, 2, 3, and 4, which contain the typical 'accretion 
states' that appeared for a given model. The characteristic 
times, are marked with arrows denoted as SI, S2 (in run 1 
and 2), S3 (only in run 2), S4 and S5 (in run 3 and 4). 

We consider the flo w to be unstable (MRI ) when the 
following relation holds (|Balbus fc Hawlevlll998t ): 

u\<-^^ (10) 
TT-^ dlnr 

where ua is the Alfven speed, Q. is the angular velocity, r 
is the distance from the BH and h is the half of the disk 
thickness. We define h as the height at which the density 
of the matter decreases by e along the vertical direction in 
comparison to the density at the equator. For run 1, 2, and 
3, the disk height can be approximated with its radius h=r 
in Eq. 1101 For run 4, we adopt the numerical value of h {h 
roughly corresponds to 1/3 r). In the vertical direction, the 
torus typically contains ~ 50, 42, 40 and 24 grid points, for 
7=5/3, 4/3, 1.2, and 1.01, respectively. In the radial direc- 
tion, there are ~ 60 and ~ 90 equaly spaced in logarithmic 
scale zones within first 20 Rs and 100 Rs, respectively. To 
find, whether the fastest growing mode of the MRI is re- 
solved with our computational mesh, we c ompute its critical 
wavelength Xcrit as defined by eq. 113 in lBalbus &i HawlevI 
(jl998h . We checked the ratio between Xcrit and Aa; (the grid 
spacing) as a function of position on the grid. If the ratio is 
larger than one, the grid resolves torus for MRI. In regions 
where the ratio is smaller than one the grid does not resolve 
the critical wavelength Xcrit- 

We have analysed the flow properties at each accretion 
state as PB03b did. For 7=5/3, we found that the accretion 
state SI corresponding to the narrow peak in Fig. [7] (panel 
for 7=5/3) appears when the high-/ torus accretes on to 
the BH. We recognize SI as 'B' state found in PB03b. Dur- 
ing SI the torus accretes because of the magneto-rotational 
instability (MRI) . We checked that the current grid does re- 
solve the fastest mode of the MRI. Xcrit /^x > 10 within 
r < 40Rs. The accretion at the peak is dominated by the 
accretion from torus, {M^/Mb reaches 0.06) and typically 
stays at this level for less than ~ 5.1 tdyn. SI appears in the 
mass accretion rate curve many times (i.e. all spikes corre- 
spond to the torus accretion). 

The accretion from torus is subsequently accompanied 
by the polar funnel inflow ('leaking') of the matter that was 
left in the polar regions after the outflow. The density in 
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the polar funnel is relatively low in comparison to the torus 
density (or in comparison to the \ow-l matter removed from 
the polar funnel by the outflow). Its value is artificially set 
by the assumed pmin parameter (see Sect. I2.2"1l . The matter 
from the polar funnel contributes to the mass accretion rate 
and sets its mean level around 0.02 (Fig. [7|. 

During state S2 in run 1, M^/M-q is very low and it 
corresponds to a state when the torus is truncated from 
the BH. The truncation is due to magnetic effects (ie., a 
'ma gnetic barrier') as in 'C state described in PB03b (see 
also lProga fc Zhandliooil and fig. 1 there). We find that the 
poloidal magnetic field is larger than the toroidal component 
for r < lO-Rs during this state. The magnetic field lines are 
stretched for r < 107?s. The MRI does not operate in the 
inner 20-Rs. The matter from the polar funnel also flows 
outward during state S2. Therefore, Ma,/ Mb drops much 
below the mean value of ~ 0.02. The duration of S2 state 
is comparable to the duration of SI and all moments with 
Mil/ Mb < 0.02 correspond to the same S2 state. 

In our run 1, we do not see an accretion state cor- 
responding to the 'D' state found in PB03b. In the 'D' 
state, Ma/Me should increase above 0.1 level, because of 
the slowly rotating matter previously pushed away by the 
coronal outflow, reentering the BH vicinity. We attribute 
lack of the 'D' state in run 1 to a relatively short duration 
of our simulations. We performed our simulation for much 
shorter time (Table [1} in comparison to PB03b (2.5 tdyn). 
In our simulations, the low-Z matter is far away (r > 500-Rs) 
from the center at the end of the simulation. 

For 7=4/3, SI state corresponds to the episode of torus 
accretion, but has a much lower amplitude than narrow 
peaks found for 7=5/3 ( SI is representative of all narrow 
peaks with amplitude of Ma/A/e ~ 0.01 in Fig. [T] panel for 
7=4/3). We find that MRI causes the torus accretion, but 
our grid resolution was not adequate to capture the fastest 
growing mode of MRI, Acrit/Ax < 1.5 within r < AQRs- 
Thus, the mass accretion rate from the torus is lower than 
it would be if the fastest mode of MRI were resolved. In 
SI state, the torus accretion is accompanied by accretion of 
the matter with density of pmin from the polar region. This 
polar funnel accretion sets the mean accretion rate at the 
level of 0.01, similar to the mean accretion rate obtained for 
7=5/3. 

During S2 state (Fig. [T] panel for 7=4/3), the torus is 
pushed away from the BH by the magnetic barrier (as in S2 
state for 7=5/3 case). 

The mentioned 'D' state appears in run 2 despite of 
a short duration of the simulation. S3 state (Fig. [71 panel 
for 7=4/3) is characterized by the the sudden increase of 
Ma/ Mb up to the level of ~ 0.2 and much longer duration 
~ SOOtdyn compared to state marked as SI. In Fig. [T] panel 
for 7=4/3, one can see only one state S3, but in Fig[l]one can 
find two broad features corresponding to state S3. During 
state S3, the low-Z matter accretes on to the BH near the 
equator. This matter has been pushed away from the BH by 
the expanding magnetic corona during the bipolar outflow 
formation but not as far as in run 1. Therefore the weakly 
rotating matter may reenter the region close to the BH in a 
much shorter time-scale. 

In the intermediate case of 7=1.2 (run 3), the mass 
accretion rate remains relatively high during the whole sim- 
ulation (Fig. [T] panel for 7=1.2) due to almost persistent 



low-/ polar radial inflow. Therefore, contrary to run 1 and 
2, it is difficult to distinguish any characteristic patterns in 
the mass accretion rate curve. S4 marks Ma/ Mb peaking 
at 0.2, when the torus accretes along with the low-Z matter. 
However, in this case the grid does not resolve the fastest 
growing mode of the MRI for r < 20i?s • The torus accretion 
is dominated by the slower MRI modes. The main contribu- 
tion to Ma comes from a polar inflow of the low-Z matter. 
The mass accretion rate is slightly higher than average be- 
cause the matter from polar regions simultaneously accretes 
from both, N and S sides (S4). State S4 appears only a few 
times during run 3. For most of the time, the matter in the 
polar regions accretes from one side only. S4 is similar to 
the S3 in run 2, but during S3 low-Z matter approaches the 
BH in the equatorial plane rather than from the poles, like 
in S4. 

In S5 state, the mass accretion rate is dominated mainly 
by the torus, because almost all matter is removed from the 
polar funnel due to the episodic coronal outburst described 
in 13.51 At that time Ma/ Mb becomes lower than 0.01. In 
run 3, the high-Z torus never undergoes the truncation be- 
cause the poloidal magnetic field is very weak and does not 
form a barrier for accretion. For 7 = 1.2 S5 appears two 
times at the end of the simulation. Note that S5 is similar 
to state SI in run 1 and 2, but here the torus accretes in- 
efficiently, therefore S5 appears as a dip (not as a spike like 
in e.g. run 1) above the mean level of accretion. 

For 7=1.01 (run 4), Ma/ Mb remains relatively high 
(Fig. [T] panel for 7=1.01) for most of the duration of the 
simulation. The low-Z matter accretes quite persistently (as 
in run 3) along the torus surface (S4), on its N or S sides, 
or both. Although MRI operates in the torus, again the ac- 
cretion fastest mode is not resolved with the simulation. 
Thus, the torus accretion rate is underestimated. In run 4, 
the torus never experiences a truncation. The low-Z matter 
inflow is only occasionally stopped by the turbulent corona 
but such episodes are short since the outflow is weak (S5). 



4 DISCUSSION AND CONCLUSIONS 

In this article we present results from the axisymmetric 2- 
D MHD simulations of a low angular momentum accretion 
flow on to the BH. Our work is an extension of MP08, 
which in turn wa s initiated by PB03a and PB03b (see also 
iProga et al]|2003l ). Here, we concentrate on investigation of 
effects of gas pressure on the flow structure and evolution. 
In particular, we study how the flow properties change as 
a function of the adiabatic index, 7. We also compare our 
MHD runs to their HD counterparts presented in MP08. 

We find that in the early phase of the evolution, the 
MHD runs behave very similar to those in MP08. Namely, 
the matter with Z > Zcr forms a torus near the equator. 
As expected, the torus thickness decreasing with decreasing 
7. The similarity in the MHD and HD initial evolution, is 
a consequence of our assumption of a initially very weak 
magnetic field. At later times, however, the MHD torii do 
not persistently accumulate matter near the BH as their HD 
counterparts do but rather accrete on the BH due to MRI 
and magnetic braking. In addition, the MHD torii are far 
more turbulent and also slightly geometrically thinner than 
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their HD counterparts. However, most importantly contrary 
to the HD torii, the MHD torii produce bipolar outflows. 

In MP08, we noticed that for 7 = 1.01, the torus is 
excited by the inflowing material and shows small amplitude 
quasi-periodic oscillations in the direction perpendicular to 
the equator. In the MHD runs, we do not observe clear torus 
oscillations in the power spectrum of the mass accretion rate 
or in any other way . This difference is due to the fact that 
in the MHD simulations the high-/ torus is more turbulent 
and does not have a well defined outer boundary what would 
separate it clearly from the inflow (such a boundary is very 
clear in the HD simulations - see the bottom row of panels 
fig. 9 in MP08). 

One of the properties of the innermost torus that is 
common for the HD and MHD simulations is the cylindrical 
distribution of the angular velocity in torus for all 7. 

The bipolar outflow and angular momentum transport 
by magnetic effects transforms the structure of the flow from 
a polar-inflow/equatorial outffow (or equatorial accumula- 
tion) in the HD models into a polar-outflow/equatorial ac- 
cretion in the MHD models. However, the torus outflow is 
strong only in models with 7=5/3, 4/3 and 1.2. For 7=1.01, 
the outflow is very weak and cannot propagate throughout 
the inflow. 

In the MHD runs with a strong torus outflow, the to- 
tal mass accretion rate of the gas on to the BH is reduced 
because Ma is dominated by the polar funnel accretion of 
the low- / gas which the torus outflow suppresses. The mass 
accretion rate is also more variable in the MHD simulations 
than in the HD ones. This variability is caused mainly by 
mixing of the low- and high-/ gas. As we showed in MP08, 
the degree of the mixing increases with the distances from 
the center for 7=4/3 and 1.2. The mixing causes episodes 
of an enhanced polar funnel accretion in the HD models 
(bursts) . In the MHD models, we observe similar bursts not 
only for 7=4/3 and 1.2 but also for 7=5/3. In both HD and 
MHD models. Ma can increase by a factor of few. For 7=1.2 
in the MHD run, we do not observe very strong flares seen 
in the HD counterpart. 

Our new simulations for 7=5/3 are consistent with the 
previous results obtained by PB03b. We found however, 
some minor differences between PB03b and our simulations. 
These differences occur during the initial phase of the evo- 
lution which we attribute to the differences in the initial 
conditions described in § [2] In addition, we do not observe 
here the episodes of enhanced accretion of the low-/ mat- 
ter (call by PB03b the accretion state D). This difference 
is due to a shorter duration of our simulations compared to 
PB03b's simulations (see the previous section). 

The torus outflow is not the only outflow we flnd in 
our simulations, i.e., we see also a highly coUimated outflow 
drive n by the centrifugal and magnetic forces. As found by 
l|Prog a 2(3Q5), this outflow appears regardless of the high-/ 
torus formation a nd is produced from the low-/ inflow (see 
IProga et al.ll200^ ). The 'outflow- from-inflow' does not sig- 
nificantly infiuence the torus evolution. If the matter in the 
polar region were not rotating (e.g. a 8 distribution of / were 
a step function at the outer boundary) this outflow would 
not form as it was the case in PB03b's simulations. 

The most important new result in the present simula- 
tions is that for a geometrically thin disc model (7 ~ 1) the 
torus outflow is strongly suppressed. Generally, the struc- 



ture of the torus and its magnetic corona depend on 7. For 
7=5/3, gas pressure is signiflcant and the flow is hot and 
the torus corona extends far from the equator. For 7=1.01, 
the magnetic corona above the disc is highly turbulent, and 
it is small and gravitationally bound. The magnetic field in 
the corona is less tangled in the fiows for 7 5/3 than the 
fiow with 7=1.01. 

As mentioned above, the torus outfiow changes the mass 
accretion rate. For 7=5/3, M^/Mb dropped suddenly about 
one order of magnitude after the outfiow develops. This drop 
in the accretion rate is due to depleting of the polar funnel. 
We observe the same behaviour of Ma/A^B in the run with 
7=4.3. For 7=1.2, the outfiow is significantly weaker and is 
often one-sided. For 7=1.01, the torus outflow is very weak 
and is stalled by the inflow. The mass accretion rate remains 
high for the most of the simulation time in comparison to 
PB03b runs and our runs for 7=5/3 and 4/3. 

In runs 1 and 2, we observe the accretion state 'C which 
occurs when the torus is truncated by the magnetic barrier. 
This state was first identified by PBQ3b. However, in the 
runs 3 and 4 (7=1.2 and 1.01) we observe the stationary 
torus accretion (called by PB03b as state B) with the state 
C never occurring. This indicates that for low 7, the role of 
the magnetic field is reduced not only in producing a torus 
outflow but also in changing the torus accretion. We note 
that the reduced role of magnetic fileds may be also due 
to resolution of our models. The MRI cannot be resolved 
with 140 X 100 points on mesh for 7 < 5/3, because the 
disk thickness decreases with 7 (half of the disk thickness 
h ~ 2/3r for 7 = 4/3 and 1.2 and h ~ l/3r for 7 = 1.01). 
We plan perform simulations with higher resolution in the 
near future. 

The direct comparison of the simulations with the ob- 
servations of BHB and AGN spectra are beyond the scope 
of this work. Therefore, we can only make some remarks on 
the observational consequences of our models. For example, 
the strong torus outfiow for 7 5/3 and the suppression of 
the outflow 7=1.01 may be relevant to the observed tran- 
sition from the hard/low state to the soft/high state. The 
episodic outburst found in models with 7=1.2 and 1.01, may 
correspond with the flares seen in a very high states in BHB. 
To relate our results with observations, we plan to carry out 
calculati ons of broad-band sp e ctra p redicted by our models 
foUowing Mos cibrodzka et al.l (|2007l ). 

Our models are axisymmetric, which may affect the re- 
sults. In axisymmetry we cannot model the toroidal fleld 
MRI and a poloidal magnetic fleld can dissipate according to 
the antidynamo theorem (Moffatt 1978). We note that even 
in the HD limit, the axisymmetr ic models differ fro m non- 
axisymmetric ones. For example. IJaniuk et al.l (|2008l ) found 
that in fully 3-D simulations the inner torus processes. It 
would be important to check if the torus will also precess 
in the MHD limit and what consequences it would have on 
the outflow formation and the polar-funnel accretion. Recent 
work on outflows driven by radiation from a precessing disc 
show that the flow geometry can be dramatically changes 
although the mass accretion rate can be similar to that in 
the non-precessing disc case yKjirosawa_& Proga 2008). We 
plan to address these issues in our future work. 

Future simulations should also include radiativ e cool - 
ing and heating processes. Recentlv. iFragile fc Meieil (|2009l ) 
pointed out that the radiative cooling may be of the same or- 
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der as artificial cooling caused by lack of the treatment of ki- 
netic and magnetic dissipative processes in non-conservative 
codes (as one used by us). Our models also do not resolve the 
geometrically thin discs. Currently, global simulations of the 
magnetized geometrically thin disc is beyond computational 
abilities. 

We acknowledge support provided by the Chandra 
awards TM7-8008X (M.M. and D.P.) and TM8-9004X 
(D.P.) issued by the Chandra X-Ray Observatory Center, 

which is operated by the Smithsonian Astrophysical Obser- 
vatory for and on behalf of NASA under contract NAS 8- 
39073. 
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Figure 1. Each panel shows the evolution of the mass accretion rate for the models in Table [Tl with different 7 parameter. 7=5/3 
(upper-left), 4/3 (upper-right), 1.2 (lower-left), and 1.01 (lower-left). The arrows N and S indicate the moments when the polar outflow 
forms in the Northern and Southern direction, respectively. The time is given in units of dynamical time at the inner boundary of the 
flow tdyn =595 s. 
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Figure 3. The total flux of total energy as a function of radius for 7=5/3 (solid), 4/3 (dotted), 1.2 (dashed), 1.01 (dot-dashed line). 
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Figure 4. The figure shows two-dimensional structures of the density and velocity fields (left panels) along with the angular momentum 
contours (middle panels), and the structure of the magnetic field parameter log(/3) (right panels), at the moments right after the 'one- 
sided polar outfiow' is formed. Each panel shows inner 200 Rs- Panels from top to bottom show outflow for different values of 7 for 
t=2800, 3400, 3450, 7500, for 7=5/3, 4/3, 1.2, and 1.01 from top to bottom, respectively. The angular momentum contours are range 
from 0.3 to 1.5 (in units of (cr) from symmetry axis toward equator with step of 0.3. Solid contours indicate I < la, and dotted contours 
I > IcT- Dotted contour of log/3 shows where log /3 changes from positive to negative value. 
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Figure 5. Same as Fig. |4]at the moments right after the polar outflow is formed on the opposite side of the torus. Panels from top to 
bottom show outflow for different values of 7 for time moments t=3500, 4000, 6000, 8300, for 7=5/3, 4/3, 1.2, and 1.01 from top to 
bottom, respectively. The value of angular momentum and log(/3) as in Fig. |4] 
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Figure 6. The evolution of density overploted with velocity field, specific angular momentum contours, and plasma parameter (3 are 
shown during an coronal outburst in run 4. The time increases from upper to bottom panels, t=9500, 10000, 10500, 11000, 13000 tdyn- 
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Figure 7. Four panels show the evolution of the accretion rate for different 7 parameter. t(jyn=595 s. 
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